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ABSTRACT 

We propose a novel method to extract the imprint of gravitational lensing by Pop III stars in the light 
curves of Gamma Ray Bursts (GRBs). Significant portions of GRBs can originate in hypernovae of Pop 
III stars and be gravitationally lensed by foreground Pop III stars or their remnants. If the lens mass is 
on the order of 10 2 — 10 3 M Q and the lens redshift is greater than 10, the time delay between two lensed 
images of a GRB is rj Is and the image separation is rj 10 /ias. Although it is difficult to resolve the two 
lensed images spatially with current facilities, the light curves of two images are superimposed with a 
delay of rs 1 s. GRB light curves usually exhibit noticeable variability, where each spike is less than Is. 
If a GRB is lensed, all spikes are superimposed with the same time delay. Hence, if the autocorrelation 
of light curve with changing time interval is calculated, it should show the resonance at the time delay 
of lensed images. Applying this autocorrelation method to GRB light curves which are archived as the 
BATSE catalogue, we demonstrate that more than half of the light curves can show the recognizable 
resonance, if they are lensed. Furthermore, in 1821 GRBs we actually find one candidate of GRB lensed 
by a Pop III star, which may be located at redshift 20 — 200. The present method is quite straightforward 
and therefore provides an effective tool to search for Pop III stars at redshift greater than 10. Using this 
method, we may find more candidates of GRBs lensed by Pop III stars in the data by the Swift satellite. 

Subject headings: gamma-ray burst — gravitational lensing — PopIII stars 



1. INTRODUCTION 

The recent observation of the cosmic microwave 
background by Wilkinson Microwave Anisotropy Probe 
( WMAP) suggests that the reionization of the universe 
took place at redshifts of 8 < z < 14 (Spergel et al. 2003; 
Kogut et al. 2003; Page et al. 2006; Spergel, D., et al. 
2006). This result implies that first generation stars (Pop 
III stars) were possibly born at z > 10 if UV photons emit- 
ted from Pop III stars are responsible for cosmic reioniza- 
tion. However, there is no direct evidence that Pop III 
stars actually formed at z > 10. Obviously, it is impos- 
sible with current or near future facilities to detect the 
emission from a Pop III star at such high redshifts (e.g., 
Mizusawa et al. 2004). But, gamma ray bursts (GRBs) 
can be detected even at z > 10, if they arise there. GRBs 
are the only currently available tool for probing first gen- 
eration objects in the universe. If one uses the data of 
absolute magnitude for GRBs with known redshifts, one 
can expect that more than half GRBs are detectable if 
they occur at z > 10 (Lamb & Reichart 2000). Roughly 
3000 GRBs have been detected to date, but redshifts have 
been measured only for 30 GRBs (Bloom et al. 2003), 
among which the most distant is GRB 000131 at z = 4.5 
(Andersen et al. 2000). But, if empirical relations be- 
tween the spectral properties and the absolute magnitude 
are used, the GRBs detected to date may include events at 
z > 10 (Lloyd-Ronning et al. 2002; Yonetoku et al. 2004; 
Murakami et al. 2005). In addition, recently a new GRB 
satellite, Swift, has been launched (Gehrels et al. 2004). 
Swift is now accumulating more data of GRBs at a high 



rate. Recently, GRB 050904 detected by Swift, in terms 
of metal absorption lines and Lyman break, turns out to 
have occurred at z = 6.295 (Kawai et al. 2006). 

The discovery of the association between GRB 030329 
and SN 2003dh has demonstrated that at least a portion 
of long bursts in GRBs are caused by collapse of massive 
stars (Kawabata et al. 2003; Price et al. 2003; Uemura 
et al. 2003). On the other hand, Pop III stars are ex- 
pected to form in a top-heavy fashion with the peak at 
1OO-1O 3 M in the initial mass function (IMF) (e.g., Naka- 
mura & Umemura 2001 and references therein). Also, the 
theoretical study by Heger & Woosley (2002) suggests that 
Pop III stars between 100M Q and 140M© may end their 
lives as GRBs accompanied by the core collapse into black 
holes. Heger et al. (2003) estimate, assuming the IMF by 
Nakamura & Umemura (2001), that 5% of Pop III stars 
can result in GRBs. In the context of cold dark matter 
cosmology, more than 10-30% of GRBs are expected to 
occur at z > 10, assuming that the redshift distributions 
of GRBs trace the cosmic star formation history (Bromm 
& Loeb 2002). Thus, observed GRBs highly probably con- 
tain GRB originating from Pop III stars at z > 10. 

The firm methods to measure redshifts are the detec- 
tion of absorption and/or emission lines of host galaxies of 
GRBs (e.g., Metzger et al. 1997), or the Lya absorption 
edge in afterglow (Andersen et al. 2000). However, these 
methods cannot be applied for all GRBs, but have been 
successful to determine redshifts only for 30 GRBs (Bloom 
et al. 2003). Instead, some empirical laws have been ap- 
plied to much more GRBs. They include a variability- 
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luminosity relation (Fenimore & Ramirez- Ruiz 2000), a 
lag-luminosity relation (Norris et al. 2000), _E p -luminosity 
relation (Amati et al. 2002), and the spectral peak energy- 
to-luminosity relation (Yonetoku et al. 2004). Apply- 
ing these relations to GRBs, the redshift distributions of 
GRBs are derived (Fenimore & Ramirez-Ruiz 2000; Nor- 
ris et al. 2000; Schaefer et al 2001; Lloyd- Ronning et al. 
2002; Yonetoku et al. 2004). Some analyses conclude that 
a portion of GRBs are located at z > 10. However, it 
is still controversial whether such an indirect technique is 
correct or not. 

In this paper, we propose a novel method to constrain 
the redshifts of GRBs that may originate from Pop III 
stars at z > 10. In the present method, the effects 
of the gravitational lensing by Pop III stars are consid- 
ered. The lensing of GRBs is considered for the first 
time by Paczynski (1986, 1987), who proposed the pos- 
sibility that a soft gamma-ray repeater is produced by 
gravitational lensing of a single burst at cosmological dis- 
tance. Also, Locb & Perna (1998) first discussed the mi- 
crolensing effect of GRB afterglows, and Garnavich; Loeb 
& Stanek (2000) found the candidate microlensed after- 
glow (GRB 000301C). The rates of such events are further 
discussed from theoretical points of view (Koopmans & 
Wambsganss 2001; Wyithe & Turner 2002; Baltz & Hui 
2005). Blaes & Webster (1992) argue the method to de- 
tect cosmological clumped dark matter by using the prob- 
ability of detectable GRB lensing. Ncmiroff et al. (1993) 
and Marani et al. (1999) search for the compact dark 
matter candidate using actual GRBs data obtained by the 
Burst and Transient Source Experiment (BATSE) satel- 
lite on the Compton Gamma Ray Observatory satellite. 
They focus on large mass lenses up to 10 6 M Q , which cause 
the delay time-scale of several tens-100 s. On the other 
hand, Williams & Wijers (1997) investigate the influence 
on GRB light curve of the millisecond gravitational lens- 
ing caused by each star in a lensing galaxy. In addition, 
Nemiroff & Marani (1998) argue that it is possible to place 
constraints on the cosmic density of dark matter, baryons, 
stars, and so on, by microlensing by stellar mass objects. 
In the present method, we focus on the gravitational lens- 
ing by Pop III stars. If the mass of Pop III stars is on the 
order of 10 2 — 1O 3 M and the redshift is greater than 10, 
the time delay between two lensed images of a GRB is « 1 
s. Quite advantageously, this time delay is longer than the 
time resolution (64 ms) of GRB light curves and shorter 
than the duration of GRB events, which is several tens 
to 100 sec for long bursts. Thus, we can see the super- 
imposed light curves of two lensed images. The present 
method seeks for the imprint of gravitational lensing by 
Pop III Stars in GRB light curves. We attempt to extract 
the imprint of lensing by calculating the autocorrelation 
of light curves. 

In this paper, we assume a standard ACDM cosmo- 
logical parameter: H a = 70 km s _1 Mpc -1 , Qm = 0.3, 
f^A = 0.7, and £lb = 0.04. The paper is organized as fol- 
lows: In §2, the formalism of gravitational lensing and the 
estimation of time delay between two images are provided. 
In §3, the method to find the evidence of lensing by Pop III 
stars is proposed. Also, we demonstrate the potentiality 
of the present method for artificially lensed GRBs, and de- 
scribe how to determine the redshifts of lensed GRBs. In 



§4, we apply this method to 1821 GRB data obtained by 
BATSE, and find a candidate of lensed GRB. §5 is devoted 
to the conclusions. 

2. GRAVITATIONAL LENSING 

We consider a GRB lensed by a foreground Pop III star. 
Here, we presuppose the lens model of a point mass. The 
Einstein ring radius gives a typical scale of gravitational 
lensing, which is expressed as 

. _ / 4GMl Dls V /2 
° E = { c 2 DosDoJ ' 
where G is the gravity constant, c is the speed of light, 
Ml is the mass of a lens object, and -Dls, -Dos, and Dol 
are respectively angular diameter distances between the 
lens and the source, the observer and the source, and the 
observer and the lens. A point mass lens produces two 
images with angular directions of 
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where (3 is the angle of lens from a line-of-sight to the 
source. We hereafter express the image with 9 > 9e by 
image 1, and that with 9 < 9e by image 2. The brightness 
of the image 1 and the image 2 are respectively magnified 
by 

A h2 = \[(i + 4/- 2 ) i /2 + (i + 4.r 2 r i /2 ± 2], (3) 

where / = (3/9e- Thus, the image 1 is brighter than the 
original one, while the image 2 is fainter. In the case of a 
lens of Pop III star, the Einstein radius is estimated as 

'-^V'V^pa. (4, 



VlO 3 M 



where D = DosDol/ D^s- Obviously, this angular sepa- 
ration is impossible to resolve by current facilities. Hence, 
we can just observe the superposition of two images. 

However, the light curves of two images are superim- 
posed with a time delay caused by the gravitational lens- 
ing, as shown in Figure 1. The arrival time of signals for 
a lensed image is expressed as 



t(9) 



(l + z L )D os D ol 
c -Dls 



-(9-[3) 2 -V(9) 



(5) 



where z L is the redshift of lens object, and \I/ is so called 
lens potential. For the point mass lens model, 'J is ex- 
pressed as 

.Dls 4GM Ll 



*(*) = 



-ln|0/6> c 



(6) 



DosDol 

where 9q is constant (Narayan & Bartelmann 1997). 
Then, the time delay between two images is given by 

At(z L , M L , /) - t(9 2 ) - t(0i) oc M L (1 + z L ). (7) 
It should be noted that At is determined solely by the mass 
and redshift of lens, regardless of the source redshift. In 
other words, the time delay places a constraint just on the 
lens, not on the source. However, if the lens redshift (zl) 
is determined, it gives the minimum value of the source 
redshift (zs) since zs must be higher than zl- 

Fig. 2 illustrates the relation between the time delay At 
and the magnification ratio between image 1 and image 2, 
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assuming the lens redshift of 50. This figure shows that 
the lens with > 10 4 M Q yields the time delay longer than 
the standard GRB duration, if the typical delay time-scale 
is assessed by / w 1. (Note that, as shown later, if / be- 
comes larger than unity, the ratio of magnification becomes 
smaller and therefore the contribution of image 2 becomes 
difficult to extract. Also, if / becomes smaller than unity, 
the probability of lensing goes down. ) On the other hand, 
the lens with < 10M© leads to At shorter than the time 
resolution of light curves, and therefore the information of 
delay is buried. The mass scale of 10 2 — 10 3 M Q expected 
for Pop III stars gives 10 _1 s < At < 1 s, which is longer 
than the time resolution and shorter than the GRB du- 
ration. Hence, this mass range appears to be suitable for 
extracting the time delay information. 

However, the actual GRB light curves generally exhibit 
variabilities with time-scale shorter than At. Thus, it is 
not straightforward to extract the time delay information. 
To demonstrate this difficulty, we show the light curve of 
GRB 930214 and the artificially lensed light curve in Fig- 
ure 3, where Ml = 10 3 M Q , zl = 50, and / = 0.5 are 
assumed. The time delay is At = 1.0 s in this case. This 
figure clearly shows that if we observe only the superim- 
posed lensed light curve, it seems impossible to recognize 
by appearance that this light curve is lensed. Hence, we 
invoke a new technique to discriminate a lensed GRB from 
unlensed one. 

3. AUTOCORRELATION METHOD 
3.1. Theory 

We pay attention to the fact that all spikes in a light 
curve are individually lensed. Then, many pairs with time 
separation of At appear in the light curve, as schemati- 
cally shown in Figure 4 (b). To detect those pairs, we em- 
ploy the autocorrelation method (e.g., Geiger & Schneider 
1996). The autocorrelation, C(6t), is defined as 
Y,iI{ti + 5t)I{ti) 



C{St) 



(8) 



where I(U) is the number of photons contained in a i-th 
bin in the GRB light curve. If there are pairs with At, 
the autocorrelation (8) is expected to show the resonance 
"bump" around At, as shown in Figure 4 (c). Then, we 
can evaluate the time delay by the existence of this bump. 

3.2. Robustness 

The autocorrelation method is simple and well defined, 
but the issue we should check is its applicability for the 
actual GRB light curves. To test the robustness of this 
method, we produce artificially lensed light curves for 
GRBs in BATSE archived data, and calculate the auto- 
correlation. We use 1821 light curves in the BA TSE cata- 
logue with the time resolution of 64 ms 4 . Unless otherwise 
specified, we adopt the data of T 90 , where Tgo is defined 
by the duration such that the cumulative photon counts 
increase from 5% to 95% of the total GRB photon counts 
(Kouveliotou et. al. 1993, Koshut et. al. 1996). Then, 
the summation in equation (8) is taken in the range of 
Tgo — St. But, if the data in Tg start with the bins whose 
time resolution is worse than 1024 ms, we neglect those 
low resolution bins. 

4 http : //cossc. gsfc.nasa.gov/batse/BATSE_Ctlg/duration. html 



In Figure 5, the resultant autocorrelation is shown for 
10 GRB light curves. In each panel, a thin solid line repre- 
sents the autocorrelation for the original light curve, while 
a thick solid line is the autocorrelation for the artificially 
lensed light curve, where Ml = 10 3 M Q , zl = 50, and 
/ = 0.5 are assumed, the same as Figure 3, and there- 
fore the time delay of lensed images is At = 1 s. We 
can see that there is no bump in C(St) for the original 
light curve, whereas a bump emerges around St = 1 s 
in C(St) for the artificially lensed light curve. Note that 
C(St) for the artificially lensed light curve is stronger than 
C(St) for all of the original light curve, owing to the am- 
plification by gravitational lensing. C(St) for the artifi- 
cially lensed light curve is fit by the polynomial of 8th 
degree. With using the best fit polynomial F(St), we 
define the dispersion, a, of the autocorrelation curve by 
a 2 = ^j = i[C(Stj) — F(5tj)] 2 /n, where n is the number of 
bins. The levels of ±3<x are shown by dashed lines. The 
zoomed view around the bump of C(St) for the artificially 
lensed light curve is also shown in each small panel. For 
these GRBs, bumps exceeding 3cr appear if lensed, cor- 
responding to the time delay between two lensed images, 
At. 

3.3. Dependence on f 

In fact, not all GRB light curves exhibit bumps in C(St) 
when lensed. The fraction of GRBs which show bumps 
exceeding 3(7 in C(St) depends on the value of / = (3/0e- 
Also, the time delay At and the magnification A\^ de- 
pend on /. To demonstrate this, we show the dependence 
on / of the autocorrelation for the artificially lensed light 
curve of GRB 930214, in Figure 6. It is clear that if / is 
larger, the bump is suppressed. This is because the con- 
tribution of image 2 becomes smaller with increasing /. 
In this GRB, the bump over 3a disappears at / > 1.0. 
The value of / at which the bump disappears differs in 
each GRB. Using 220 GRB data, which are used in Fen- 
imorc & Ramirez- Ruiz (2000), we obtain the fraction of 
GRBs which show bumps over 3er in C{5t) as a function 
of /. In this analysis, At — 1.0 s is assumed. The resul- 
tant fraction is shown in Figure 7 (a). For / = 0.5, about 
a half of GRBs show bumps exceeding 3er in C(5t). In 
contrast, the cross-section of the gravitational lensing is 
proportional to / 2 , which is also shown in Figure 7 (a). 
We evaluate the probability density of GRBs exhibiting 
bumps over 3cr by multiplying the fraction for which a 
bump appears by f 2 . The normalized probability density 
against / is shown in Figure 7 (b). As a result, the proba- 
bility is peaked around / = 1 and the standard deviation 
corresponds to A/ 0.25. It is noted that this probabil- 
ity density is found to be hardly dependent on the value 
of At. 

3.4. Optical depth 

Here, we estimate the optical depth of gravitational lens- 
ing by Pop III stars. If we assume that the fraction a of 
baryonic matter composes Pop III stars at z > zni, then 
the optical depth is given by 



r z s c dt 
r(z s )= n L (z L )o- L —dz L , 

J z m 



'dz L 



(9) 
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where <7l is the cross-section given as ox = ^(-Dol^e) 2 
and til(zl) is the number density of the lens objects given 
as 

3H r^b /1 x3 



(10) 



Then, the optical depth (9) is calculated as 



r(z s ) = -afl b 
5 



(1 + z s ) 5/2 + 1 



(1 
(1- 



Zs )5/2 _ 
)5/2 



1 



(1 + Z„l) 5 /2 _ 1 



ln(l + z s ) 
ln(l + z m )| (11) 



in the Einstcin-de Sitter universe (Turner, Ostriker, & 
Gott 1984; Turner & Umemura 1997). Here, we assume 
a = 0.1 and zm = 10. The resultant optical depth of 
Pop III star lensing is shown in Fig. 8. Since t(zs) is 
the probability that a source is located inside the Einstein 
ring radius (/ < 1), it is not the probability of bump de- 
tection. The probability of bump appearance, p(f), is a 
decreasing function of /, as shown in Fig. 7 (a). Since the 
optical depth that / is in the range of [/, / + df] is given by 
dr(zs) — r(zs)df 2 = r(zs)2fdf, the probability of bump 
detection is given by 

P(z s ) = r 2fp(f)r(z s )df. (12) 
Jo 

The resultant bump detection probability is also shown in 
Fig. 8. From this figure, the probability turns out to be 
« 0.001 for zs = 20 — 40. If we take into account that 
more than 10-30% of GRBs occur at z > 10 (Bromm & 
Loeb 2002), the expectation number of bump detection 
for lensed GRBs is assessed to be one in a few thousand 
GRBs. 



4. A CANDIDATE FOR GRB LENSED AT Z 



60 



4.1. Data analysis 



As shown above, the autocorrelation of intrinsic light 
curves exhibits no bumps for almost all GRBs. But, a few 
in 1000 GRBs might show bumps in C(5t) even for intrin- 
sic light curves. Hence, we calculate the autocorrelation of 
all GRB light curves available in BATSE catalogue, which 
amount to 1821 GRBs. As a result, we have found one 
candidate, GRB 940919 (BATSE trigger number 3174), 
in which a 3a bump in C(St) appears. The light curve 
and the autocorrelation of this GRB is shown in Figure 
9. As seen in panel (b), a bump exceeding 3a appears at 
At = 0.96 s. 

4.2. Statistical significance 

To check the statistical significance of a bump in C(6t) 
of GRB 940919, we make a test with mock light curves. 
Here, we generate mock light curves using a smoothed cor- 
relation function that dose not show any bump, and inves- 
tigate whether bumps appear in correlation functions just 
from pure statistical fluctuations. 

From Wiener-Kihntchinc theorem, the power spectrum 
of light curves are given by 

|/MI 2 = ^E^(Wi + ^)], (is) 

where T denotes the Fourier transformation and uo = 
2n/8t. We can generate mock light curves by the inverse 



Fourier transformation of I(lu). Here, in order to add fluc- 
tuation to I(ijS), we take random Gaussian distributions, 
where \I(oj)\ is the standard deviation and the phase is 
random in the range of [0, 2n]. Then, the mock light curve 
is given by 

I(t) = = £ w |/»| cos(-«/> - cot), (14) 

where J(w) is a random sample in the Gaussian distribu- 
tion, and <f> is the random phase shift from to 2tt. We 
produce 2000 mock light curves using a correlation func- 
tion, and recalculate the autocorrelation C(8t) by equation 
(8). As a result, we have found that no bump higher than 
3a appears in the C{8t) of 2000 mock light curves. A part 
of the recalculated C(8t) are shown in Figure 10. Thus, it 
is unlikely that a bump in the correlation arises as a result 
of pure statistical fluctuations. 



4.3. Light curve decomposition 

As a further test for the lensing of GRB 940919 light 
curve, we attempt to decompose the light curve, assuming 
that it is the superposition of two lensed light curves with 
At = 0.96 s, and analyze the decomposed light curves. 
The decomposition is made by the following recurrence 
formula; 



itot(i) = h{t) + h{t), 



h(t) = ^h(t 



At), 



(15) 



(16) 



where 7 t ot (t) is the observed intensity, and I\ (t) and I2 (t) 
are intensities for image 1 and 2, respectively. If these two 
equations are combined, I\(t) can be expressed by 



N f A V 

j=0 ^ 1 / 



(17) 



The summation is taken in the range of t—jAt ^ 0, where 
t = is the starting point of Tg . Then, we can derive also 
h{t) by equation (15). In Figure 11, the decomposed light 
curves are shown in the case of / = 1 which is the most 
probable case as shown in §3.3. The application of this 
decomposition method for the finite amount of data does 
not guarantee that the light curve is successfully decom- 
posed into two lensed light curves. Hence, to check the 
validity of this decomposition method, we calculate the 
cross-correlation of two decomposed light curves by 



C c (8t) 



J2ih(U)l2(U + St) 



(18) 



where 8t is the time shift. The result is shown in Figure 
12. As clearly shown in this figure, the cross-correlation 
is peaked when St accords with At = 0.96 s. Also, each 
decomposed light curve shows no bump higher than 3a in 
the autocorrelation for a reasonable range of /. Hence, we 
can conclude that the light curve of GRB 940919 is suc- 
cessfully decomposed and is likely to be the superposition 
of two lensed light curves. 
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4.4. Redshift estimation 

Here, we constrain the redshift of the lens object. As 
shown in equation (7), At just determines Ml(1 + Zl), ex- 
cept for /. If the probability density against / (Fig. 7 (b) 
) is applied, we can derive a suitable range for Ml(1 + zl)- 
The dark gray region in Figure 13 represents the suitable 
range for At = 0.96 s. On the other hand, the hatched 
region shows the mass range of Pop III stars obtained by 
Nakamura & Umemura (2001). Combining these two re- 
gions, the allowed redshift of the lens object is at least 
zl ~ 10. If / = 1 is adopted, the redshift ranges from 
zl « 20 to zl ~ 200, where the most probable one is 
z L m 60. 

Nonetheless, there still exists another possibility that 
the lens object is as massive as ~ 10 4 M Q located at zl — 0, 
as seen in Fig. 13. Loeb (1993) and Umemura et al. (1993) 
suggested that relic massive black holes are candidates for 
such an object. Sasaki & Umemura (1996) place a con- 
straint on Hbh from the UV background intensity and the 
Gunn-Peterson effect in the context of a cold dark mat- 
ter cosmology. They find that the black hole mass density 
might be as low as ClBn/Qb ^ 10~ 3 . Therefore, the ex- 
pected number of bump detection for massive black holes 
is by two orders of magnitude lower than that for Pop III 
stars. 

5. CONCLUSIONS 

To place constraints on the redshifts of GRBs which 
originate from Pop III stars at z > 10, we have proposed 
a novel method based on the gravitational lensing effects. 
If the lens is Pop III stars with 10 2 - 1O 3 M at z > 10, 
the time delay between two lensed images of a GRB is « 1 
s. This time delay is longer than the time resolution (64 



ms) of GRB light curves and shorter than the duration of 
GRB events. Therefore, if a GRB is lensed, we observe 
the superposition of two lensed light curves. We have con- 
sidered the autocorrelation method to extract the imprint 
of gravitational lensing by Pop III stars in the GRB light 
curves. Using BATSE data, we have derived the probabil- 
ity of the resonance bump in the autocorrelation function, 
which is an indicator for the gravitational lensing. Apply- 
ing this autocorrelation method to GRB light curves in the 
BATSE catalogue, we have demonstrated that more than 
half of the light curves can show resonance bumps, if they 
are lensed. Furthermore, in 1821 GRB light curves, we 
have found one candidate of GRB lensed by a Pop III star 
at z « 60. The present method is quite straightforward 
and therefore provides an effective tool to search for Pop 
III stars at redshift greater than 10. Although the num- 
ber of GRBs with available data is 1821 in this paper, the 
Swift satellite is now accumulating more GRB data. If the 
present method is applied for those data, more candidates 
for GRBs lensed at z > 10 may be found in the future. 
These can provide a firm evidence of massive Pop III stars 
born at high redshifts. 
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Fig. 1. — Schematic diagram of lensed GRB light curve. The dotted line represents the original light curve, while the dashed and dot-dashed 
lines represent the light curves of image 1 and image 2, respectively. The superimposed light curve of lensed images is shown by the solid line. 
The intensities are normalized by the maximum intensity of original light curve. 
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Fig. 2. — Relation between the time delay and the magnification ratio (A2/A1) between two images. The redshift of lens object is fixed to 
50. The curves show lens masses of 10-10 5 A/q with an increment of 1 order of magnitude. In the right vertical axis, / = 0/6&, is also shown 
corresponding to the magnification ratio. The time resolution of BATSE data, At = 64 ms, is represented by the leftmost gray region. 
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Fig. 3. — The original light curve of GRB 930214 obtained by BATSE and artificially lensed light curves are presented. The lens is assumed 
to be a Pop III star with Ml = 1O 3 M0 at zl = 50, while the source is located at zg = 51. Lensed image 1, image 2, and the total light curves 
are shown. The time delay between two images is At = 1.0 s under the assumption of / = 0.5. 




Fig. 4. — (a) Schematic figure of a spiky GRB light curve, (b) The light curve as the superposition of two lensed light curves with the time 
delay of At. (c) Autocorrelation function calculated for light curve (b). 
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Fig. 5. — Autocorrelation against St for 10 GRB light curves. In each panel, a thin curve shows the autocorrelation for the original light 
curve, while a thick curve is that for the artificially lcnsed light curve. Here, M L = 1O 3 M , z L = 50, and / = 0.5 are assumed and therefore 
the time delay of lensed images is At = 1 s. A dotted curve is the best fitting for the autocorrelation for lensed light curve, and two dashed 
curves show ±3cr level from the best fitting. In each panel, a zoomed view around a correlation bump is also shown. 
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Fig. 6.- 
panel (b) 
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- Dependence on / = /3/6e °f the autocorrelation for the artificially lensed light curve of GRB 930214, which is the same GRB as 
in Fig. 5. The time delay is At = 0.1 s. The meanings of curves are the same as Fig. 5. It is seen that the bump is weaker for 




Fig. 7. — (a) A dot-dashed curve shows the probability of exhibiting bumps exceeding 3cr. Here, At = 1.0 s is assumed. A dashed line 
is f 2 , which is proportional to the cross-section of gravitational lens, (b) The probability density of exhibiting bumps against /. The gray 
scales represent lcr, 2<j, and 3cr around / = 1 (from dark to light). 



Imprint of Gravitational Lensing 



11 



0.006 





12 



Hirose et al. 



1.04 



1.02 




0.5 1 1.5 2 2.5 3 



8t[s] 



Fig. 10. — Autocorrelation for mock light curves. A thick curve is the autocorrelation for the original light curve (GRB 000421). Other 
thin dashed curves are a part of autocorrelations for 2000 mock light curves based on the same GRB. 
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Fig. 13. — The dark gray region shows the suitable range of M-l(1 + zl) for GRB 940919. Here, 3<r range of / (see Fig. 7b) is adopted. 
The central thick solid line corresponds to the case of / = 1, which is most probable. The hatched region is the mass range of Pop III stars 
by Nakamura & Umemura (2001). 



